Symmetry breaking by the sea of Dirac-Landau levels in graphene 
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The quantum Hall states of graphene have a filled sea of lower Dirac-Landau levels. The short 
ranged SU{4) symmetry breaking interactions can induce a staggered polarization of the sea of 
Dirac-Landau levels. We study this effect in the extended Hubbard model on a honeycomb lattice 
using mean field variational wavefunctions. We find a valley symmetry broken, anti-ferromagnetic 
spin ordered phase at v — ±1 when the on-site interaction is dominant. Our mean field solution is 
consistent with the recently reported experimental results of Z. Jiang et. al.[l] 



The low energy physics of graphene is described by 
quasiparticles satisfying the massless Dirac equation [2] 
[3]. There are two species of Dirac fermions per spin 
making a total of four species of quasiparticles. In the 
non-interacting limit, these four species are degenerate 
resulting in an SU(4) symmetry of the non- interacting 
theory. The dominant interaction is the long-range part 
of the Coulomb interaction which is also SU(4) symmet- 
ric. 

The quantum Hall effect has been observed in graphene 
[1] [2] [3] [4] [5] . At low magnetic fields (~ 1-10 T), 
plateaus are observed at {h/e 2 )<JH = 4(n + 1/2) con- 
sistent with there being four degenerate Dirac-Landau 
levels. At higher fields (~ 20-45 T), this degeneracy 
is lifted [5]. It is completely lifted in the n = level 
resulting in extra plateaus at (h/e 2 )<JH = —1,0,1 and 
partially lifted for n = ±1 yielding extra plateaus at 
(h/e 2 )(TH = ±4. Tilted field experiments [5] indicate 
that the plateau at (h/e 2 )aH = ±4 is due to the Zeeman 
splitting of the Landau levels. Recent experiments by 
Z. Jiang et. al. [1] show that the transport gap at the 
{h/e 2 )<JH = ±1 (y = ±1) plateau does not depend on the 
parallel component of the magnetic field and is propor- 
tional to y/B±. These experiments hence indicate that 
the v = ±1 plateau corresponds to a state with zero net 
spin polarisation. These plateaus are attributed to the 
lifting of the sub-lattice degeneracy of the n = Landau 
level. 

It has been shown that the long range part of the 
Coulomb interaction causes the SU(A) symmetry to be 
broken spontaneously [6] [7] [8] [9]. The resulting gap is 
proportional to y/B±. The exact pattern of the symme- 
try breaking is not determined by the ST/ (4) symmet- 
ric long range part of Coulomb interaction but by other 
symmetry breaking interactions [10] . These could be the 



Zeeman term, short range lattice scale interactions [11] 
[8] [12] and disorder [13]. The effects of the short range 
lattice scale interactions can be studied in the extended 
Hubbard model on the honeycomb lattice in the pres- 
ence of a magnetic field. The previous studies of this 
model [8] [11] find a spin polarised state at v = ±1. How- 
ever, this is not consistent with the experiments discussed 
above [1]. 

In this work, we are mainly concerned with the physics 
of v = ±1 states and we show that the sea of Dirac- 
Landau levels plays a significant role. We write down 
mean field trial wavefunctions that incorporate a stag- 
gered SU (4) polarization of the n/0 Landau levels and 
compute their energies analytically in a systematic a/l c 
expansion [9], where a (= 2.45 ^4°) is the lattice spac- 
ing and l c (f=a 40 A° at 45 T) is the magnetic length. 
We then find that when the on-site repulsion is domi- 
nant, the ground state has anti- ferromagnetic ordering 
with broken valley symmetry. It has zero net spin po- 
larisation and hence is a strong candidate for the state 
found in the experiments [1]. 

The hamiltonian we consider is, 
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hi = ^2 a c\ a c i( j is the number operator and </>y is the 
phase due to presence of the magnetic field. The con- 
tinuum model describing the low energy physics can be 
derived in the standard way. The continuum effective 
hamiltonian is, 



H = Jd 2 x ¥(v F a.U)^ - ^ ^(*V a *) 2 + (*r 2 cr a *) 2 + i(*t a y a a ^) 2 ^j + ((*t$)2 _ (|r r ^) 



(2) 



Where II is the covariant derivative. The kinetic en- ergy term has an SU (4) internal symmetry whereas the 
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interaction terms break it down to [Z^ x U(l)) va u ey <E) 
SU(2) sp i n . The remanent of the SU(2) valley symmetry 
of the non-interacting theory corresponds to 
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where, r is the Dirac (sub-lattice) index, 77 the val- 
ley index and a the spin index. Note that the Z2 does 
not commute with the U(l). The Z2 part of the valley 
symmetry corresponds to the interchanging of the valley 
indices. The U (1) part corresponds to the conservation of 
the number of quasiparticlcs belonging to the two valleys 
separately. This is a symmetry only of the low energy 
sector as scattering from one valley to the other is a high 
momentum transfer process. There are also no Umklapp 
processes which contribute to this process. 

Our variational wavcfunctions generalise those written 
down by Yang et. al. [7] to include the polarisation of the 
Dirac sea. They are constructed from the eigenfunctions 
of the following massive single particle Dirac equation, 
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where m is a real number and Q a matrix constructed 
from a basis of four orthonormal SU(A) spinors, 
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for the case where k of the four n 
The variational states are, 
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$"'(x)x p being the eight component spinor eigenfunc- 
tions of the hamiltonian in Eq.(5). The cutoff on the 
number of Landau levels, N c , is obtained by matching the 
total number of states with that of the lattice. We have, 
N c = (27r/V3)(^/ a 2 ). Note that in the limit m -> 0, 
our variational states reduce to those of Yang et. al. [7]. 
The ground state manifold is U{A)/{U(k) x [7(4— k) for 
all to. The single particle energy levels are shifted by the 
presence of the mass term as shown in Fig.(l). 

For the computation of the expectation value of the 
hamiltonian, we need to compute the coincident two 
point correlation function 
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FIG. 1: (a) shows the energy levels for non-interacting Dirac 
particles in magnetic field for graphene.(b) shows the shift in 
one particle levels due to mass and the filling of one particle 
states for v = — 1. 



Assuming that to ~ (a/Z c ) 2 we put to = rh(a 2 / \2-kID) . 
This assumption will be justified later when we solve for 
to. r n can be computed to leading order in a/l c to be, 
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Equation (10) shows that the mass term induces a stag- 
gered SU(4) polarization of the n 7^ Landau levels. 
Since the Dirac index is the same as the sub-lattice in- 
dex, the order parameters for staggered 51/ (4) order are 
(4>A a *) = Tr /3A a r, where A a are the SU(A) generators. 

Using the two point coincident correlation function we 
can compute the expectation values of hamiltonian in 
terms of Q the matrix. The kinetic energy density, E t , 
interaction energy densities Ey and Ejj arc computed to 
be, (here we have retained the terms upto leading order 
in a/l c and dropped the terms that are independent of 
to and Q) 
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where, n = 1 + (2m /i) and i = (t/2)(3 v / 37r) 1 / 2 . Note 
that the m — > limit is obtained by putting fx = 1 i.e. 
the contribution of filled sea of Dirac-Landau levels is 
neglected to the mean field energy. 

We now specialise to the case when v = — 1 (k = 1) 
where a single n = Dirac-Landau level is filled. The 
ground state manifold is parametrised by a single SU (4) 
spinor which specifies the SU (4) polarization. The v = 
1 (k = 3) case is related to the former by a particle-hole 
transformation. We use an explicit paramcterisation, 
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This corresponds to a linear superposition of an elec- 
tron with valley index +, spin polarization h\ and valley 
index — , spin polarization — n®. 7 and Q specify the 
relative amplitude and phase of the superposition. 

When the 7 is either or ir the many body state has 
a definite number of electrons in each valley. The n = 
level electrons are localised in one sub-lattice with arbi- 
trary spin, corresponding to charge and spin ordering. 
These are the states discussed in previous work [8] [11]. 
They have the U(l) part of the valley symmetry unbro- 
ken and the Z2 part broken. 

When < 7 < 7T, the many body state does not have 
a definite number of electrons in each valley. The total 
number however, remains a good quantum number. Fur- 
ther, if 7 7^ ±7r/2, then the average number of electrons 
in each valley is not the same and there is charge order- 
ing. In this case the full valley symmetry, Z 2 x U(l), is 
broken. When 7 = ±7r/2, the average number of elec- 
trons in the two valleys are equal, the Z2 symmetry is 
unbroken and the U(l) is broken. 

In general the state will have a non-zero total spin 
polarization except when 7 = ±7r/2 and hi — fi2- In this 
case there is anti-ferromagnetic spin order and no charge 
order. 

We now evaluate the energy density in terms of our pa- 
rameters by substituting Eq.(14) in equations (6,12 and 
13). The total mean field energy density (E) is the sum 
of kinetic (E t ), nearest neighbour interaction (Ey), Hub- 
bard interaction (Ejj) and Zeeman term (Ez)- 



previously [8] [11]. It is interesting to note that when the 
Dirac sea contributions arc neglected (fi = 1), Ey is in- 
dependent of the Q. This was also noticed by Alicca 
and Fisher [8] . They found that making the interactions 
slightly non-local picks out the 7 = 0, it state. 

The Hubbard term, Eu, does not contribute to the 
mean field energy if we negelect the contributions from 
the sea of filled Dirac-Landau levels. This can be quali- 
tatively understood in the n — subspace. Since the two 
valley species live on the the two distinct sub-lattices, if 
only one n = level is occupied, we can have any SU (4) 
polarization without any double occupancies. However, 
when the n/0 levels are taken into account, Ejj is min- 
imized when 8 = and 7 = ±7r/2. The ground state 
is thus anti-ferromagnetic spin ordering (SDW) with the 
U(l) valley symmetry broken. 

We need to minimise this energy density Eq.(15) with 
respect to the variational parameters 8, 7 and m. Let us 
first consider the case when we ignore the Zeeman energy. 
The energy density minimizes for 7 = 0, tt (independent 
of 8) and 7 = ±ir/2 (8 = 0). The minimum energy is 
obtained at, 
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7=|(-|) if 3V-U<0; 8 = 



Thus at large V we have the Z2 broken, U(l) unbroken, 
spin polarised state and at large U we have the Z2 unbro- 
ken, £7(1) broken, anti-ferromagnetic phase. The transi- 
tion between them is first order and occurs at 3V = U (at 
mean field level). The phase diagram is shown in inset 
of Fig. (2) 

The minimization of energy density Eq.(15), yields 
same two phases in the V-U space discussed above. The 
masses in the two phases are given by 
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The transition remains first order and the line separat- 
ing two phases is given by the solution of 
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g is the Zeeman parameter which is 60K for B ~ 45T 
and cos# = hi ■ n-i 

Ey is minimized at 7 = or 7 = 7r. Thus the nearest 
neighbour interaction picks out the Z2 broken, f (1) un- 
broken spin and charge ordered state (CDW) discussed 
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Note that to leading order in a/l c , the phase boundary is 
independent of the magnetic field. In general we expect 
a weak magnetic field dependence. 

The resulting phase diagram for the v = ±1 is shown 
in Fig. (2). Thus at large U the Dirac sea polarization 
effects can drive the system into a (7(1) valley symmetry 
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broken anti-fcrromagnetic phase. Since the total spin 
polarization of this state is zero, the gap will have no 
dependence on the parallel component of the magnetic 
field in this phase, which is consistent with the recent 
experiments [1]. We therefore have a possible mechanism 
to explain the tilted field experiments at filling factors 
v = ±1 [1]. It has been argued earlier, the Hubbard U 
drives the system towards anti-ferromagnetic order [14]. 
It can be seen from the Fig(2) that the estimate for the 
critical Hubbard J7(~ 15 eV) is large. This number could 
change by including higher order corrections in a/l c and 
the fluctuations about the mean field. Correlations may 
also be important even at the integer fillings [15] [16]. 
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Thus we have shown that if among the short-range 
part of the interactions, the Hubbard interaction is dom- 
inant, then graphene has a U(l) valley symmetry broken 
phase with zero spin polarization at v = 1. The gap in 
this phase is proportional to y/B±. This phase is there- 
fore consistent with the experimental results reported [1] 
where it is seen that the gap varies as y/B± and is inde- 
pendent of -By. 

There will be three gaplcss collective modes in this 
phase. Two corresponding to the ST/ (2) — ► U(l) spin 
symmetry breaking and one corresponding to the U(l) 
valley symmetry breaking. The topological defects in 
this phase are clearly interesting objects to study. 

We have also analysed the v = phase in de- 
tail. The state is now paramctcriscd by two orthogonal 
SU(A) spinors with ten parameters corresponding to the 
U(4)/(U{2) x U{2)) ground state manifold. The results 
will be reported in a longer forthcoming publication along 
with more details of the results reported in this letter. 

Acknowledgement: We thank G. Baskaran for many 
useful discussions. 



FIG. 2: The figure shows the phase diagram for v = ±1 in 
the presence of the Zeeman term. The inset shows the phase 
diagram for the case when the Zeeman term is neglected. The 
CDW region, Z2 part of the valley symmetry is broken. In 
the SDW region, (7(1) part of the valley symmetry is broken. 
All energies are in eV. 

As we mentioned earlier, the dominant interaction in 
graphene is the SU (4) symmetric long range part of the 
Coulomb repulsion. As shown by Yang et. al. [7], this 
interaction favours ST/ (4) "ferromagnetic" order by the 
direct exchange mechanism. The SU (4) ferromagnetic 
order is quantified by the value of (W. ^ $ a i) . As can be 
seen from equation (9), this is non-zero for the variational 
states we are considering. It is in fact independent of m 
because it gets contributions only from the n = levels. 
Thus the dominant gap in graphene will come from the 
the Coulomb exchange interaction which is proportional 
to y/BZ ■ 

The variational states also have an SU (4) "anti- 
ferromagnetic" order corresponding to a staggering at 
the scale of the lattice spacing. This is quantified by the 
value of (^ V a^ifa')- Equation (10) shows that this is the 
component that is enhanced by a non-zero value of to. 
However at the scale of l c over which the exchange mech- 
anism operates, the effects of this order which varies at 
the scale of a will be small. The dominant gap will there- 
fore be the "ferromagnetic" gap discussed above which is 
proportional to s/B±. 
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